Reentrant charge ordering caused by polaron formation 
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Based on a two-dimensional extended Hubbard model with electron-phonon interaction, we have 
studied the effect of polaron formation on the charge ordering (CO) transition. It is found that for 
fully ferromagnetically ordered spins the CO state may go through a process of appearance, collapse 
and reappearance with decreasing temperature. This is entirely due to a temperature-dependent 
polaron bandwidth. On the other hand, when a paramagnetic spin state is considered, only a simple 
reentrant behavior of the CO transition is found, which is only partly due to polaron effect. This 
model is proposed as an explanation of the observed reentrant behavior of the CO transition in the 
layered manganite LaSr2Mn2 07. 

PACS numbers: 71.30.-Fh, 71.38.-Fi, 75.50Cc 



Electron crystallization was first proposed by Wigner, 
who considered an electron system in a uniform posi- 
tive background. It was pointed out that at sufficiently 
low densities the Coulomb energy of electrons dominates 
their kinetic energy and then the electrons form a lat- 
tice. Such a Wigner lattice is experimentally realized in 
semiconductor heterostructures [||. 

Wigner's scenario is also important for solids Elec- 
tron or hole crystallization, i.e., charge ordering (CO) is 
extensively observed in real materials, even at high den- 
sities. The Verwey transition in magnetite Fe304 is a 
typical example of a real space ordering of Fe^+/Fe'^+ 
species Furthermore hole ordering is found in rare- 
earth pnictides like Yb4As3 because of intrinsic small hy- 
bridization between Yb-4f and As-4p orbitals The 
compound a'-NaV2 05 which was originally proposed as 
another inorganic spin-Peierls material is now known to 
have a CO of V4+/V5+ at 34K §]. CO also appears in 
a variety of colossal magnetoresistance manganites, e.g., 
in Lai_2;Caa;Mn03 for x > 0.5 In all the above ex- 
amples, only one CO transition is observed. 

Recently, a melting of the CO state on decreasing the 
temperature, i.e., reentrant behavior, has been found in 
layered mariganite LaSr2Mn207 by electron and x-ray 
diffraction It was proposed that a CO state ex- 

ists characterized by superlattice reflections appearing 
at ~ 210K, but collapsing again below lOOK. Subse- 
quently, more precise measurements have found that sur- 
prisingly the superlattice intensities start growing again 
below about 50K meaning the reappearance of the CO 
state at lower temperatures That is to say, in the 
whole temperature region the CO state goes through a 
process of appearance, disappearance and reappearance. 

It is the purpose of this letter to propose a mecha- 
nism which can explain this intricate behaviour. It is 
known that the essence of CO transition is competi- 
tion of the inter-site Coulomb energy and kinetic energy 
of electrons. Therefore a tempting idea is to invoke a 
temperature-dependent electron bandwidth such that the 



kinetic energy is small compared to the Coulomb energy 
at some temperatures (favoring CO state), but dominates 
at other temperatures (favoring homogeneous state). 
Such a mechanism can be realized by including electron- 
phonon interaction, i.e., considering polaron formation 
since the effective hopping of a polaron is drastically 
affected by the number of thermally excited phonons. 
Indeed there is a lot of evidence for polaron formation 
in manganite materials, for example in Lai.2Sri.8Mn207 
|lO| . Very recently, small polaron ordering was di- 
rectly observed in the CO phases in La(Ca)Mn03 and 
Pr(Ca)Mn03 0. Theoretically the CO transition has 
already been studied within the polaron context [ p^ , but 
more profound physics has been ignored. This will be ad- 
dressed in the following. Rather than directly discussing 
material aspects we focus on a general interesting prob- 
lem: how is the CO transition itself affected by a simulta- 
neous polaron formation? The inherent evolution of the 
ordered polaron state with temperature is presented and 
rich reentrant behavior is found. 

Our model is described by the following Hamiltonian 
in a two-dimensional square lattice: 

H = -t'^ (cj ,^Ci+5,cr + h.c.) + U'^ Ui^Uii 

i,5,(T i 
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where ^{ci^a) denotes creation (annihilation) opera- 
tor for an electron at site i with spin a, Ui is defined 
as Ui = n^i -t- nil with ~ cl^Ci,^, t is nearest- 
neighbouring hopping parameter and U, Vi , V2 represent 
on-site, n.n. and n.n.n. Coulomb repulsion, respectively. 
5 — x.y are unit vectors along the x and y direction and 
T] = X ± y. The first four terms describe an extended 
Hubbard model. Here we consider the electron-phonon 
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interaction in the following way. We assume that be- 
tween arbitrary two n.n. sites there exists another kind 
of ion, e.g., oxygen, which vibrates along the bond di- 
rection. Such a vibration will be coupled to the electron 
density difference of two neighboring sites, as shown in 
the above fifth term [Q. Then m is the oxygen ion 
mass, uj is vibration frequency, Ui^s represents the dis- 
placement of the oxygen ion located between sites i and 
i + S. The vibrations have been assumed dispersionless 
for simplicity. Our electron-phonon coupling is not of 
an usual on-site form which may be used to model the 
coupling with the apical oxygen vibrations since experi- 
ment on layered manganite Lai.2Sri.8Mn207 showed ||ic| ] 
that polaronic distortions are constrained within the per- 
ovskite layers . For further discussion, we rewrite the 
electron-phonon interaction in the boson representation: 
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where g = f / \j2muj. We set ?i = fcs = 1 throughout and 
all energies are in units of K. 

The total number of sites is N and for the electron 
band we assume quarter-filling. We will be interested 
in CO with wavevector (1/2, 1/2) and consider different 
occupancies on the two sublattices A and B of a bipartite 
lattice. First we treat the electron-phonon interaction 
with a combined Lang-Firsov (LF) transformation \J\ — 
exp[- Yl,i,s al^^i.s - - "i)] 111; ^nd squeezing 

transformation U2 = exp[7 1]i,5&<5^i,<5 ^ '^i*^^ 7 

(> 0) a variational parameter |16[ . Then averaging on 
equilibrium phonon states p7| , we obtain the following 
effective electronic Hamiltonian: 

iieSS = ^(r + 1/t)7V/2 - i ^(4,c,+5,, + h.c) 

i.S, a 
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where r — exp(— 47) and i = t exp[— 5ar coth(ijj/2r)] 
with definition of dimensionless coupling constant a — 
juP'. It is obvious that electron-phonon interaction 
leads to two effects. One is that the bare hopping t 
is reduced to t. The narrowing effect which is largely 
overestimated by the LF method alone in adiabatic limit 
(ujjt << 1) is weakened by the squeezing effect. 
Our emphasis here is the fact that i itself increases with 
decreasing temperature because the phonon occupation 
decreases until only the contribution coming from zero- 
point vibrations is left. Such a temperature-dependent 
renormalizcd hopping is very important for reentrance 
of CO transition, as will be seen later. Another effect is 
phonon-mediated interaction between electrons which fa- 
vors CO state, as described by the last term above. In our 



above treatment, the equilibrium positions of vibrating 
oxygens have been displaced, e.g., (ui^s) (n^ — rii+s), 
which means the oxygens will be repelled to unoccupied 
sites once CO forms. In the following we concentrate on 
the CO transition for Hamiltonian (||). It is complicated 
if we simultaneously consider charge and spin degrees of 
freedom. We will therefore treat only cases with explic- 
itly given spin configuration. 

First, we consider the fully ferromagnetical (FM) spin 
ordered case which is equivalent to the case of spinless 
electrons at half filling for which the Hamiltonian (||) , up 
to a constant, is simplified as 



Heff = c^(t + l/T)iV/2 - f^(a|c,+5 + h.c.) 

(4) 
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with Ui = cjci. In mean-field (MF) treatment and 
with the assumption {m) = l/2±x,i G A, B {x rep- 
resents deviation from homogeneous electron distribu- 
tion), the above Hamiltonian can be diagonalized with 

energy bands ±ej- ^j: — + (AVx)'^ with tj: — 

4<cosA;i cosfc2 (-vr/2 < ki, k2 < 7r/2) and V = V12 + 
2auj {V12 — Vi — V2 > 0) in the CO phase. It is easy to 
find that in this case the chemical potential /i = for all 
temperatures, see also Ref. jl^]. Note that it is only the 
difference V12 of intersite Coulomb energies that enter 
the quasiparticle energy. Since Vi, V2 are well screened 
and of the order 0.3eV in manganites ||2^ the difference 
V12 as estimated from neighbour distances is of the or- 
der lO'^K. For this reason the CO transition may occur 
in an observable temperature range. More precisely it 
is the difference of Madelung energies (with long range 
Coulomb interactions included) between the charge ho- 
mogeneous and the CO phases that sets the scale of the 
transition temperature. Using the density of states for 
the spectrum s^: 
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where K is the complete elliptic function of the first kind 
and z = \J — (AVx)'^ j^t we finally get the following 
self-consistent equation for order parameter x 
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It can be easily proven that this equation always has a 
solution with nonzero x at zero temperature no matter 
what are the values of parameters t, a, V12 . In the current 
half-filled spinless case this means that the system must 
take on CO state at zero temperature. We will refer to 
this fact later on. For given initial values for t,Vi2,uJ 
(for example, t = 4000,^12 = 400, w = 200), we give 
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the solution of the above equation in the phase diagram 
T vs a, shown in thick Hne of Fig. [|. The parameter 
T (or 7) is determined by minimization of ground state 
energy of the homogeneous state. Its minor variation 
with parameter a is ignored since it does not change our 
quahtative results. An average value r = 0.12 is adopted. 




FIG. 1. Phase diagram T vs a for the case of fully FM spin 
order with parameters t = 4000, V12 = 400, u = 200. The 
thin line is a hypothetical CO phase boundary for a constant 
t — texp( — 5ar) (see text). The inset shows the tempera- 
ture-dependent order parameter x. The dotted line, solid line, 
dashed line correspond to a = 0.60, 0.65, 0.68, respectively. 
The multiple reentrant behavior is obvious for a = 0.65. 

From the phase diagram it can be clearly seen that 
there is an intermediate region for a where the system 
goes through a process of homogeneous state, CO state, 
homogeneous state and CO state with decreasing tem- 
perature. The corresponding variation of the order pa- 
rameter X with temperature is shown in the inset of Fig. 
|l| by a solid line. The physical explanation is as follows. 
At high temperature, the hopping is drastically reduced 
so that CO state is favored to appear at a relatively high 
temperature (It is usually quantitatively overestimated 
due to our simple treatment) . On decreasing temperature 
the effective hopping increases. And most importantly, 
it increases so fast that its corresponding CO transition 
temperature decreases even faster than the temperature 
itself. Therefore the CO state cannot be reached any 
more and the system returns to the homogeneous state. 
However, at even lower temperature, the system is sure 
to recover to CO state since the ground state of the sys- 
tem is charge ordered, as refered above. Actually such a 
conclusion is strongly favored by other methods beyond 
MF theory Physically one has a half-filled band 

of spinless electrons with perfect nesting which has a CO 
instability for arbitrary small inter-site Coulomb interac- 
tion at zero temperature. In this sense it is inevitable 
that the system finally enters the CO state. We want to 
emphasize, the multiple reentrant behavior found here is 
entirely due to a temperature-dependent polaron band- 
width. For comparison, a hypothetical phase diagram 
without the temperature-dependent factor in hopping t, 
i.e., i — iexp(— 5ar) is shown as thin line in Fig. |^. 
Obviously no rcentrance of CO happens now. 



The phase diagram with its reentrant behavior is sim- 
ilar for different values of parameters. However it is ap- 
propriate to discuss the influence of the optical phonon 
frequency to. It is obvious that the renormalized hop- 
ping t is nearly unchanged with the temperature in the 
region about T < lu/2 since the phonons are difficult to 
excite. This means that the temperature-dependent ef- 
fective hopping which we emphasized here is mainly ob- 
served in the region T > uj/2. Therefore, roughly speak- 
ing, if we increase the frequency the reentrant region will 
as a whole move to higher temperatures. It should be also 
clearly pointed out that our temperature dependence of 
the effective hopping t is not quantitatively exact. How- 
ever the precise form of i{T) is not crucial for the above 
reentrance behavior as has been checked with a func- 
tional dependence of t{T) less pronounced than the ex- 
ponential behavior given before. 




FIG. 2. Same to Fig. |l| for the case of PM spin order with 
parameters t = 4000, Vo = -2000, = 200. In the inset the 
solid and dashed lines correspond to a = 3.7, 4.1, respectively. 
Only one reentrant transition is obvious for a = 3.7. 

Now, we discuss the case including spin and assume 
a paramagnetic (PM) spin state with = (nii) = 

{ni)/2. For moderate U and with similar steps as before 
we can obtain self-consistent MF equations for the order 
parameter x and chemical potential /i: 
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F^{z)K{\/l - z2) dz , (6) 
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with functions F±{z) = l/{l + exp[-{4:^{izy + {Vxy + 
fi)/T]) ± 1/(1 + exp[(4V(iz)2 + (Vx)'^ - /i)/r]) and the 
new definition V = Vo + 3aw (Vb = V12 - U/8). At 
zero temperature it is easy to obtain the condition: 
V/t > const « 1.2, under which a solution with nonzero 
X exists. Obviously a must exceeds a threshold to reach 
CO state as long as initially Vq < 1.2t. This is a major 
difference from the fully FM spin ordered case. The nu- 
merical solution for Eqs. (3) and (^ is summarized in 
the phase diagram of Fig. 2 with a thick line (we have 
set the same value for t as before). Due to different zero 
temperature properties now only possibility of one reen- 
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trance is found as shown by the sohd hne in the inset. For 
comparision, we also give the hypothetical phase diagram 
with i = texp(— 5ar) by a thin line in Fig. |[ Reentrant 
behavior is also found. Actually, our hypothetical case is 
qualitatively equivalent to considering a T vs V12 (with 
fixed U) phase diagram when electron-phonon interac- 
tion is not included. Even without electron-phonon in- 
teraction, reentrant behavior is also possible in some V12 
region in the PM spin case, see also Ref. |^^. However, 
the quantitative reentrant behavior is much different with 
and without polaron effect. 

It is natural to think about a crossover between the 
two extreme spin cases, PM and complete FM, i.e. to 
consider the case of general incomplete FM spin order. 
Since the fully FM spin order can be considered as a 
limit of on-site repulsion U —^ 00, such a crossover has to 
be treated beyond MF which will be investigated in the 
future. We predict that even in the intermediate case the 
multiple reentrant behavior is controlled by the polaron 
effect. The extended Hubbard model itself is expected to 
show at most one reentrant transition. This conjecture 
is obtained from an interpolation between Fig.l and 2, 
see the thick and thin phase boundaries respectively. 

Finally we give a possible explanation of the exper- 
imental observations in manganite LaSr2Mn207 in the 
context of our theory. This material is composed of 
Mn02 bilayers and there is one Cg electron every two 
manganese ions. When double-exchange mechanism is 
considered, the initial hopping t should be substituted 
by t{cos{9/2)) [6: angle between two t2g spins on n.n. 
sites), which depends on local spin order. There is strong 
evidence that in this material in-plane FM spin or- 
der has already been well developed within the Mn02 
layers at room temperature, before it exhibits true long 
range order. This is also concluded from the observation 
of FM in-plane spin waves even at twice the 3D tran- 
sition temperature So a good assumption is con- 
sidering that the spins (for both t2g and eg electrons) 
are perfectly aligned in the whole temperature region as 
far as a single layer is concerned. For the same reason 
the dependence of t on the temperature resulting from 
double-exchange mechanism is ignored. Then we can use 
our theory with fully FM spin order to explain the CO 
scenario of LaSr2Mn207. The reentrant behavior seen in 
Fig. is qualitatively the same as found experimentally 
. Moreover the experimental observation of reappear- 
ing superlattice intensities at lower temperatures found 
in Ref. |^ is naturally explained within our theory. We 
point out, the multiple reentrant behavior is also possible 
from our theory even if the assumption about perfect FM 
spin order is relaxed. Of course, full explanation on the 
experiments needs to consider other factors, especially an 
orbital degree of freedom ||2^. However, we think that 
polaron formation plays an essential role for the observa- 
tion of reentrant CO transitions in this material. 

In conclusion, in a 2D extended Hubbard model we 



have found interesting reentrant behavior of CO tran- 
sition when including the polaron effect. Due to a 
temperature-dependent polaron bandwidth the CO state 
may go through a process of appearance, collapse and 
reappearance on decreasing the temperature in case of 
fully FM spin order. When a PM spin order is considered, 
once reentrance of CO transition may happens, which is 
partly due to polaron effect. The crossover between these 
two cases is complicated and still has to be investigated. 
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